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This paper presents an aeroelastic finite-element formulation for axially loaded aerodynamic structures. 

The presence of axial loading causes the bending and torsional sitffnesses to change. For aircraft with axially 
loaded structures such as the truss-braced wing aircraft, the aeroelastic behaviors of such structures are non- 
linear and depend on the aerodynamic loading exerted on these structures. Under axial strain, a tensile force is 
created which can influence the stiffness of the overall aircraft structure. This tension stiffening is a geometric 
nonlinear effect that needs to be captured in aeroelastic analyses to better understand the behaviors of these 
types of aircraft structures. A frequency analysis of a rotating blade structure is performed to demonstrate the 
analytical method. A flutter analysis of a truss-braced wing aircraft is performed to analyze the effect of geo- 
metric nonlinear effect of tension stiffening on the flutter speed. The results show that the geometric nonlinear 
tension stiffening effect can have a significant impact on the flutter speed prediction. In general, increased 
wing loading results in an increase in the flutter speed. The study illustrates the importance of accounting for 
the geometric nonlinear tension stiffening effect in analyzing the truss-braced wing aircraft. 

I. Introduction 

Future aircraft design trends are employing high-aspect ratio wing technology to improve aerodynamic efficiency 
for reduced fuel burn. As a wing aspect ratio increases, the need for maintaining sufficient structural load carrying 
capacity is becoming more important. Cantilever wing design can only accommodate up to a certain aspect ratio 
beyond which the wing root bending moment can become too large that will impose structural limitations on a wing 
design. At the same time, wing deflection can also become too excessive that could cause aerodynamic penalty. 
Truss-braced wing aircraft concepts provide a structural solution to high aspect ratio wing aircraft designs. The long 
slender wing would employ structural bracing via the use of axially loaded strut members to provide intermediate span 
supports in addition to the wing root attachment. These struts generally support a portion of the spanload carried by the 
wing and are loaded in tension. Under a negative-g flight condition such as during a dive, a load reversal could occur 
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that could put the struts in compression. The compressive loading would require design considerations for buckling 
strength. 

Under aerodynamic loading, an axially loaded member also experiences the normal bending and torsion generated 
by aerodynamic lift force and pitching moment. Aeroelasticity of an axially loaded structure undergoing transverse 
bending can be significantly different from that with transverse bending alone. 

The SUGAR Truss-Braced Wing (TBW) aircraft concept is a Boeing-developed N+3 aircraft configuration funded 
by NASA ARMD Fixed Wing project. 1 ' 2 The TBW aircraft concept is designed to be aerodynamically efficient by 
employing a high-aspect ratio wing design. The aspect ratio of the TBW is in the order of 14 which is significantly 
greater than those of conventional aircraft wings. As a result, intermediate structural supports are required. The main 
wings are braced at approximated at mid-span by two main struts. In addition, two jury struts; one on each wing, 
provide additional reinforcement. Figure 1 is an illustration of the TBW aircraft. 



Fig. 1 - Boeing SUGAR Truss-Braced Wing (TBW) Aircraft Concept 

The main struts and jury struts can impact aerodynamics of the main wings since they create aerodynamic inter- 
ference drag that can offset the aerodynamic benefit of high-aspect ratio wing design. Thus, the TBW is designed to 
strive a balance between aerodynamic efficiency and structural efficiency. 

In this paper, we will develop an aeroelastic analysis for axially loaded aerodynamic structures such as the TBW 
main struts. 3 The presence of axial loading causes the bending and torsional sitffnesses to change. A tensile loading 
will result in an increase in bending and torsional stiffnesses. This tension stiffening is a geometric nonlinear effect 
which should be captured in aeroelastic analyses of aircraft structures that carry significant axial loads. A finite- 
element formulation is developed to model the axially loaded structures in the presence of the bending-torsion motion 
of a typical aircraft wing. A flutter analysis is performed to analyze the geometric nonlinear effect of tension stiffening 
on the flutter speed. 
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II. Elastic Analysis 


In the subsequent analysis, the combined motion of the left wing is considered. The motion of the right wing is 
a mirror image of that of the left wing for symmetric flight. The wing has a varying pre-twist angle y(x) commonly 
designed in many aircraft. Typically, the wing pre -twist angle varies from being nose-up at the wing root to nose-down 
at the wing tip. The nose-down pre-twist at the wing tip is designed to delay stall onsets. This is called a wash-out 
twist distribution. Under aerodynamic forces and moments, the aeroelastic deflections of a wing introduce stresses 
and strains into the wing structure. The internal structure of a wing typically comprises a complex arrangement of load 
carrying spars and wing boxes. Nonetheless, the elastic behavior of a wing can be captured by the use of equivalent 
stiffness properties. These properties can be derived from structural certification testing that yields information about 
wing deflections as a function of loading. For high aspect ratio wings, an equivalent beam approach can be used to 
analyze aeroelastic deflections with good accuracy . 4 

Consider an airfoil section on the left wing as shown in Fig. 2 undergoing bending and torsional deflections. Let 
(x, y, z) be the undeformed coordinates of point Q on a wing airfoil section in the reference frame D. Let po = xdi be 
a position vector along the elastic axis. Then, point Q is defined by a position vector p = Po + q where q = yd 2 + zd 3 
defines point Q in the y — z plane from the elastic axis. Then the undeformed local airfoil coordinates of point Q are 
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where 77 and q are local airfoil coordinates, and y is the wing section pre -twist angle, positive nose-down . 5 
Differentiating y and z with respect to x gives 
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Fig. 2 - Left Wing Reference Frame of Wing in Combined Bending-Torsion 

Let © be a torsional twist angle about the x-axis, positive nose-down. Let W and V be flapwise and chordwise 
bending deflections of point Q, respectively. Let U be the axial displacement of point Q. Then, the displacement and 
rotation vectors due to the elastic deformation can be expressed as 

r = t/di+Vd 2 +Wd 3 (3) 
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(j> =@di-W x d 2 + V x d 3 


(4) 


where the subscripts x and t denote the partial derivatives of 0, W, and V. 

Let (xi,yi,zi) be the deformed coordinates of point Q on the airfoil in the reference frame D and pi = v i d i + 
y \&2 +zid 3 be its position vector. Then the coordinates (x \ . y i . z i j are computed as 

pi=p + r + 0xq (5) 


where 
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Differentiating x\,y\, and zi with respect to x yields 


V la- 


1 + u x - 7 V™ + zyVx - zW xx - yyW x 

71 a- 

= 

~zy +V x -z® x -yy® 

. ZM . 


77 + w x + 7 ® a - zy © 


Neglecting the transverse shear effect, the longitudinal strain is computed as 6 
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The axial force and moments acting on a wing are then obtained as 6 
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where A P x and AMi x ^ z \ are the nonlinear terms; E is the Young’s modulus; G is the shear modulus; A is the tensile 
area; e y and e z are the offsets of the centroid from the elastic axis; I yy , I yz , and I zz are the section area moments of inertia 
about the flapwise axis; J is the torsional constant; and B\, B 2 , and are the bending-torsion coupling constants which 
are defined as 
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The slope of the twist angle y can play a significant role in structures with large twists such as turbomachinery 
blades. For aircraft wings, this effect is negligible and therefore can be neglected. Also, in general, the offset of the 
centroid from the elastic axis in the ^-direction, e z , is small compared to that in the y — direction e y . Therefore, it may 
be neglected. Then, the moments acting a wing are reduced to 
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Note that the signs of the moments are defined in the positive deflection sense such that 


M = M x di - M y d 2 + M,d 3 


(17) 
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The nonlinear terms are given by 
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where 


Consider two cases: 
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1. If the nonlinear strain component due to axial displacement (xi^ 
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If, in addition, the nonlinear strain component due to chordwise deflection 



0, then the nonlinear 
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terms are simplified to 
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III. Aeroelastic Angle of Attack 


The relative velocity of the air approaching a wing section includes the contribution from the wing elastic deflection 
that results in changes in the local angle of attack. Since aerodynamic forces and moments are dependent on the local 
angle of attack, the wing aeroelastic deflections will generate additional elastic forces and moments. The local angle 
of attack depends on the relative approaching air velocity as well as the rotation angle 0 from Eq. (4). The relative 
air velocity in turn also depends on the deflection-induced velocity. The velocity at point Q due to the aircraft velocity 
and angular velocity in the reference frame D is computed as 
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and ( p,q 7 r ) are aircraft angular velocity components in the roll, pitch, and yaw axes, and (x a ,y a ,Z a ) is the coordinate 
of point Q in the aircraft body-fixed reference frame B relative to the aircraft C.G. (center of gravity) such that x a is 
positive when point Q is aft of the aircraft CG, y a is positive when point Q is toward the left wing from the aircraft 
C.G., and z a is positive when point Q is above the aircraft C.G. 

Consider a trim problem when p = 0, p = q = r = 0. Then 
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The local velocity at point Q due to aircraft rigid-body dynamics and aeroelastic deflections in the reference frame 
D is obtained as 7 

dr (90 
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where 
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In order to compute the aeroelastic forces and moments, the velocity must be transformed from the reference frame 
D to the airfoil local coordinate reference frame defined by (jl . t] . J as shown in Figure 2. Then the transformation 
can be performed using successive rotation matrix multiplication operations as 
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for small deflections. 

The local aeroelastic angle of attack on the airfoil section due to the velocity components Vjj and in the reference 
frame D, as shown in Fig. 8, is computed as 
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Assuming that z ~ 0 and neglecting any third-order terms or higher, the local aeroelastic angle of attack is evaluated 
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as 
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Let n sa VL, w ~ VL a, and vv,- = VLa,, where a, is the induced angle of attack along the wing span. Then the 
nonlinear partial derivatives of the local aeroelastic angle of attack are evaluated as 
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We note that the time derivative of the axial displacement contributes to the aeroelastic angle of attack as a damping 
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source. Otherwise, the axial displacement U is weakly coupled to the wing aeroelasticity. 
The nonlinear angle of attack is then expressed as 


a c (x,y) = a r (x) + a e (x,y) 


where a, (x) is the rigid angle of attack and a e (x.y) is the elastic angle of attack 
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The elastic angle of attack can be linearized as a function of the static deflection as 
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Finally, if we ignore all the deflection-dependent terms, then the linearized elastic angle of attack becomes 
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IV. Aerodynamic Forces and Moments 

The wing section unsteady lift coefficient is given by 
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where k = is the reduced frequency parameter, ft) is the frequency of wing oscillations, c is the section chord in the 
streamwise direction, c/ a is the section lift curve slope, a ac is the effective elastic angle of attack at the three-quarter- 
chord point, 5 and C (k) is the Theodorsen’s complex-valued function. 8 

The section pitching moment coefficient about the aircraft pitch axis is evaluated as 


ft r /-t / t \ i ftw ft CXacC 
Cm — C,n ac + ~CL a + C ( k ) a ac \ — 

c c 2Voo 


(56) 


where c,„ ac is the section pitching moment coefficient at the aerodynamic center, positive nose up, e is the offset of the 
aerodynamic center forward of the elastic axis, and e u is the offset of the three-quarter chord point aft of the elastic 
axis. 

The section drag coefficient is expressed in a parabolic drag polar form as 


cd — c Do + K ft/. 


( 57 ) 
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where cd 0 is the section parasitic drag coefficient and K is a parabolic drag parameter. 
The force and moment coefficients can be written as 
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where c^ r , c mr , and cu r are the rigid lift, pitching moment, and drag coefficient contributions. 

The lift force, drag force, and pitching moment in the aircraft coordinate reference frame along the elastic axis are 
expressed as 
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The axial and normal force coefficients are defined as 
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A. Propulsive Forces and Moments 

In addition, the propulsive effects of the aircraft engines must be accounted for in the analysis. Both the engine mass 
and thrust can contribute to the aeroelasticity. 7 The propulsive force and moment vector are expressed in the reference 
frame D as 


( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 

(71) 


(64) 

(65) 


f = 8 (x-x e ) (F * di +F y e d 2 + F z e d 3 ) 

= 8 (x — x e ) [(— T sinAcosT — m e gsinr) di — T cos Ad? + ( T sinAsinT — m e gcosr)d 3 ] (72) 

nT = r f xf f = (x e d , -y e d 2 -z e d 3 ) x f,, (73) 
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where T is the engine thrust, m e is the engine mass, (x e ,y e ,z e ) is the coordinate of the engine thrust center such that 
y e > 0 forward of the elastic axis and z e > 0 below the elastic axis, and 8 (x—x e ) is the Dirac delta function such that 


j 8(x-x e )f(x)dx = f{x e ) (74) 

Transforming into the local coordinate reference frame and neglecting nonlinear contributions, the propulsive 
forces and moments are given by 

ft = 5 (x-x e ) (f; + F,y x + F y e V x ) (75) 

ft = 8 (x-x e ) [f; + F* (© + y) - F x e V x ] (76) 

ft = 8 (x - x e ) [Ff - F; (© + 7 ) - F x e W x \ (77) 

m x = S(x- x e ) [FyZe - F*y e + (F x y e + F y x e )W x - {F x z e + F z x e ) V v ] (78) 

m e y = 8(x- Xe) \—F x z e - F z x e + ( F x y e + F y x e ) (© + y) - (F y z e - F z y e ) V x ] (79) 

m e z = 8 (x - x e ) [F x y e + F y x e + ( F x z e + F z x e ) (© + /)- ( F y z e - F z y e ) W x \ (80) 


B. Inertial Forces and Moments 

The inertial forces and moments are due to acceleration experienced by the wing, the air volume surrounding it which 
gives rise to the apparent mass acting at mid-chord, as well as the engine mass and half of the mass of the fuselage and 
tail empennage. The acceleration results in the inertial forces and moments as follows: 


n 2 ... , 1 


f x = —mg sin r — mU tt — — cos AU tt — 8 (x) -mfgsinT 

ft = —mV t t - yPcoC 2 cosAV), + 8 (x-x e ) (-m e z e &„ - m e V„ ) + 8 (x) ( ftn f Zf& r t - ^m/V, 


/l' O /(, 

ft = —mg cosT + me cg &tt + — p TO c cos Ae,„0 ff — mW tt — —pooC COsAW tt + 8 (x — x e ) ( m e y e &it — m e W t t) 


(81) 

(82) 

(83) 


+ 5 (x) ^m/gcos T - ^ m /)’/&„ - ^m f W„^J 

JC y y 7 V A TC 'J 

m x = mge cg cosT — mftO u — —p^c e;,, cosA@ tt — —— PooC cos A® tt + me C gW tt + —pooC~e m cosAW tt 

4 128 4 

+ 8 (x - x e ) (-m e ft®„ + m e y e W,t - m e z e V„) + 8 (x) (^-^mjgyj cos V - gft® u - ^m f yfW„ + \m f ZfV,^j (84) 


1 ] 

m'y = 8 (x) ( —-mfgZfSmT + 


m‘ z = —mge cg sin V + 8 (x) ( ftn f gy f sinT- -l{ z V xtt 


(85) 

( 86 ) 


where e m is the offset of the mid-chord aft of the elastic axis; e cg is the offset of the center of mass forward of the 
elastic axis; (0 ,yy,Zf) is the coordinate of the aircraft CG such that yy > 0 aft of elastic axis and Zf > 0 above the 
elastic axis; nif is the mass of the fuselage and tail empennage; and I ft, I ft, and //. are the inertias about the left wing 
reference frame of the fuselage and tail empennage. 

The aircraft CG is located in the aircraft body-fixed reference frame B at a distance x/ ; aft of the elastic axis and zb 
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above the elastic axis at x = 0. Then by coordinate transformation into the reference frame D 


yf = jc'bcosA 

Z/ = — XflsinAsinr + zbcosT 


(87) 

( 88 ) 


The roll, pitch, yaw, and roll-yaw inertias of the fuselage and the tail empennage at the elastic axis at x = 0 are 
computed as 


// =// 
1 XX — l xx 


Iyy ~ 4 + m f ( x b + Zb) 

4 = 4 + m f x B 

4 = 4 + m fXBZB 


(89) 

(90) 

(91) 

(92) 


-T -f -f -f 

where I xx , Iy Y , I zz , and I x/ are the roll, pitch, yaw, and roll-yaw inertias at the aircraft CG. 
The inertias can be transformed into the wing reference frame D as 


iL = 


— sinAcosT — cosAcosT — sinT 


1 

X 

0 

1 

N 


— sinAcosT 

0 

jf 

l YY 

0 


— cos A cos r 

i 

X^> 

N 

0 

hz 


-sinT 


= it x sin“ A cos z r + Iy Y cos" A cos z T + I zz sin z T — It 7 sin A sin 2r 


jf 


f 


I f = 

yy 


if = 

zz 


— cos A sin A 0 


sin A sin r cos A sin T — cosT 


X 

0 

N 


0 

if 

Lyy 

0 


N 

0 

fzz 



— cos A 
sin A 
0 


= I xx cos 2 A + Iy Y sin 2 A 


1 

4 

0 

1 

N 

s a,* 


sin A sin r 

0 

if 

lyy 

0 


cos A sin r 

1 

N 

0 

Izz 


— cosT 


= lf 


XX ^ 


siir Asin z r + ll Y cos z A sin z F + I zz cos A T + I X7 sin A sin 2r 


■f 


(93) 

(94) 


(95) 


V. Aeroelastic Equations for Combined Extensional, Bending, and Torsion Motion 

A. Equilibrium Conditions 

Consider the equilibrium of forces and moments acting on a wing section. The resulting equilibrium equations are 
given by 6 



dp x 

,+A = ° 

ox 

(96) 


* I'A 3 
+ 

II 

o 

(97) 


dPy 

(98) 

dM x 

dx 

-P y W x + P z V x + m x = 0 

(99) 
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( 100 ) 


These equations become 


dM v 

-j^-PxW x +P z -m y = 0 
dM- 

— P X V X + P y + m z — 0 


dM z \ (dM v \ 

dx - P,V X + m z P X W X - my J V x + m x = 0 


d 2 M v d (P X W X ) dm v 

—^~2~ ~ j --r-~fz = 0 (103) 

dx L dx dx 

d 2 M- d (P X V X ) dm 7 

+ = 0 (104) 

dx~ dx dx 

Let U (x,t) = U (x) + u (x,t), © (x,t) = © (x) + 0 (x,t), W (x,t) = W (x) +w(x,t), and V (x,t) = V (x) + v (x,t) 

where the bar symbols denote the static deflection. The static tensile force acting on the wing section is defined as 
T = EAU X . Then, considering the linear deflection terms, we have 

P x EAu x - EAe y v xx 

M x (GJ + Tk 2 )G x 

= ' ’ (105) 

My E IyyW XX E lyyV XX 

Af - EAc y u x EI yz w xx T EI zz v xx 


d ( EAu x EAeyVxx) 


d [( Gj+Tk 2 ) e x 

dx 


= ~m x - — - - P X V X + m z )w x - - P x v x - P X V X + m z W x 


( dM y - . _\ 

+ ( - P X W X - m y J v x - 


J LV lx; - - \ - 

- P x w x - P X W X - m y j PV X 


d 2 {EIyyW XX EIy Z V XX ) 3 17ly 3 (P X W X + P X W X ) / 1 f\Q 

T~2 — Jz i 5 I 5 

dx A dx dx 

d ( EAe y u x Ely 7 w xx -\ EI zz v xx ) dm z | d ( P x v x + P x V x ) /1An . 

dx 2 ^ 

where I xx = Ak 2 . 

Since V x and V xx are generally small, so they can be neglected. Then, P x = T so that the equations above become 


d[(Gj+Tk 2 )e x 

dx 


\d(-Tey-EI yz W xx ) 1 

lx- — 7 — 1 I -rn z M 

dx 

d ( EAc y u x EI yz yv xx A E I zz y xx ) 
dx 


1 - [ d(EIyyW XX ) - 1 

-Tv x + m z W x + f XX1 -TW x -m y v x (110) 
dx 


d~ (EIyyW xx EIy Z v xx ) dm y d 7'vv x T {EAu x EAe y v xx ) W x ] 

~dx 2 ~ h+ ~d7 + dx 
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( 112 ) 


d~ ( EAeyUx EIy 7 w xx 4 EI zz v xx ) d m 7 d (T v x ) 

dx 2 y dx dx 

Finally, if the vertical bending deflection is not too large so that the nonlinear effect is not dominant, then in the 
simplest form these equations are expressed as 


d [(GJ + Tk 2 ) 0 X 

^ = -m x 

dx 

d 2 ( EIyyW xx - EI yz v xx ) _ dirty d (T w x ) 

dx 1 ' dx dx 

d~ ( EAe y u x EI yz w xx + E I 77 y xx ) dm z d ( Tv x ) 

dx 2 ' y dx dx 


(113) 

(114) 

(115) 


The aeroelastic equations of motion for axially loaded structures with no deflection-dependent terms are given by 

d ( EAu x — EAe y v xx ) 


dx 


[rn + ^p„oC 2 cos u n - 8 (x-x e ) (F^w x + F y v x ) 


(116) 


d \ (GJ + Tk 2 ) 6 X 1 ?? / , n 7 9 n , \ 

— s ~ — c m qooC~cos~ A cosT + mr k + — p«,c e“ cosA+ — -p^c cos A G tt 

dx \ 4 12s / 

- (me cg + 7 ^p 0 cC 2 e m cos w tt -8(x-x e ) [(F x y e +F y x e )w x - (F x z e + F z x e )v x ] 

+ 8 (x-x e ) (m e r 2 G tt - m e y e w„ +m e z e Vtt) + 8 (x) + \4 x G tt - \m f ZfV t , 


(117) 


d (ElyyW xx EIy z V XX ) A \ ^ / 2 ' \ \ ^ 2 , \ 

•- 9 = cl^ooCCOsA+ — [c m q«,c sinAcosA) + \ me cg + — p^c e m cosA J G tt 

- (tn + ^c 2 cos A ^ w tt + 8 (x - x e ) (-F y 0 - F x w x ) 

+ { 5 (x - x e ) [( F x y e + FyX e ) 6 - (. F y z e - F z y e ) v x ] } + 8 (x - x e ) ( m e y e G t , - m e w„) 


,1 1 \ d 

+ o (x) ( --m f yf6t, - -m f w tt J + 


s (■ x ) ( 2 J yy Wxtt 


d(T \ i 


dx 


(118) 


d ( EAc 7 u x E Iy 7 w xx + E I 77 v xx ) 

d? 


7 8 t 77 , / 7T a \ 

b< 7 ocCCOs A — (c m qooC ~ cos AsinT) — ^m+ — pooC cos AJ v u 


+ 8 (x - x e ) (Ff 6 - F x v x ) - J- { 8 (x - x e ) \(F x z e + F z x e ) G - ( F y z e - F z y e ) w x ] } 


+ 8 (x Xe) (- MeZeOtt - m e v tt ) + 8 (x) ( \m f ZfG n - \m f v t,'] - 


8 ( x )[~2 I (z v xtt 


d (Tv x 
dx 


(119) 
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B. Discretization 


The partial differential equations can be solved by weak-form solutions such as the finite-element method (FEM). 
Using interpolation functions, the displacements can be expressed as 


n 

u(x,t) = Y N u( x ) u i (0 
i=l 

o(x,t) = Y Ns ( x ) e '(0 
1=1 
n 

w(x,t) = Y,N w (x)<l)i(t) 
i=i 
n 

v(x,t) = YjNv WfiW 
1=1 


( 120 ) 

( 121 ) 

( 122 ) 

(123) 


where N u (x), Ng (x), N w (x), and N v (x) are the interpolation functions, 0,- 
n is the number of elements. 

The discretization results in a spring-mass-damper system in the form 


- 

T 


Wi w x j 

. Wi = 

Vi V x>i 


and 


n 


n n 


Y, M i X i + Y. C i*i + Y K ' X ‘ 

i= 1 i= 1 i= 1 


= o 


(124) 


where Mj is the elemental mass matrix, C, is the elemental damping matrix, K, is the elemental stiffness matrix, and 

lT 


0, vv,- 


Wr 


Vi Vv 


is the nodal displacement vector. 


The mass matrix due to the structural mass, aerodynamic mass, engine mass, and fuselage mass is given by 


Mi = M) + M \ f -(- Mj + m\ 


(125) 
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where 


Mf = 


NjN u 

0 

0 

0 


NjqNe 

N w e cgN o 

0 


N J 6 C gN w 


N 1 N 

iy w iy w 


Mf = 


f 1 ' 71 2 A 

/ — PooC"cosA 

Jo 


0 

0 

0 


7 t 7 

w 

0 

0 


0 

0 

0 

nJn v 


dx 


N e ( e ?K + 32)^0 

h 1 0 


A^g 

Nj,N w 


0 

0 

0 

ivTm, 


dx 


+ / qo°ccos z A < 



0 

0 

0 

0 


0 

^0 ccosr ^0 


N eCcosr^N v 


0 



-NJ,^X^N V 


_ 0 


-K T ^ 



+ 


O 

0 

0 

O 


0 

0 

0 

0 


0 

N'JctanA^Ne 

N'Jc tanA f^w 

ctanA^M, 


O 

-N' v T csmr^N e 

-iv' T csinr|^iV w 

-iV' T csinr^A v _ 

> 


> t/x 


Mf = m„ 


' = 2 


0 

0 

0 

_ 0 

0 

0 

0 


A^JV 0 

-Afovty, 

NJz e No 


Aq y<>N w 
N^N w 


Njz e N v 

0 

A^JM. 


NjlLNe&i 
Nj,m f yfNe 
[ 0 —NjnifZfNe 


0 

Ng m fyfN w 
Nj, m f N w +N'J l{yN w 

0 


0 

— A^g nifZfN v 
0 

<m/iV v +^ T 7^ v 


*=0 


The damping matrix due to structural damping and aerodynamic damping is given by 

C,'=C, s + Cf 
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(128) 


(129) 


(130) 
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where 


Cf = 2M/4>C^ _1 


Cf = / g^ccos 2 A < 


0 

0 

0 


Njccosr-^Ng 


— Af T 1 g£riy _mT 1 

JV w COR A df). y O 1Y w cnsA 


v w cos A dfy 


AjccosT (§*“#, 
( 

( 


w cos A 

-Nj 


££i A7 
dw, Vw 


i£a-AJ 

dw, w 


dCm 

dw x , 

dc L 

dwri 


A,. 


A, 


dc a 


dwxt w 


A, 


A/T r rn c I | ^C,,, /V J_ df m a/ 

jv 0 ccosr ^ 5 Vf ^Vvi- 5vja rv v 
A 
( 


—A, 


T 1 


w cos A 

-Nj 


(^v 


^A 


ac„ ^ , a^ »/ 
av, /v '’ + av A / v v 


o 

0 

0 

O 


0 

0 

0 

0 


0 

A' T ctanA^A 0 

A7ctanA(^A w+ |^A;) 

A7ctanA(^A v + f^A;,) 


O 

-A' T csinr^A 0 

-A; T csinr(^A w +J^A;) 

-A; T csinr(^A l ,+ ^A:) _ 

y 


(131) 


Jx 


(132) 


where £ = diag ( Ci , § 2 , • • • , Cn) is the diagonal viscous damping ratio matrix, <f> is the structural dynamic eigenvector 
matrix, and El = diag i(0 \ . o>>. . . . . 0J„ ) is a diagonal matrix of the structural frequencies. 

The stiffness matrix due to structural stiffness including the nonlinear tension-induced stiffness, aerodynamic 
stiffness, and engine thrust-induced stiffness is given by 


Kj = Kf + Kf + Kf 


(133) 


where 


Kf = 


A JEAN U 
0 
0 

-A l T EAe y N' u 


0 

A 0 T (GJ + Tk 2 )N e 
0 
0 


0 

0 

N^EIyyK + N 'J TN' w 

-N" J EI y -X 


Kf = / ^ooCCOs 2 A< 


0 
0 
0 
0 

0 

ATf = ‘ ° 
0 

0 



0 

0 

0 

0 



0 

NjccosF^fNg 

N e ccosr ^F x N 'w 

NqC COST XjN 



0 


-N J 1 I^a' 

w cosA dw x w 

-N t 1 & N 

w COS A dv x a 



_ 0 

-A V T d XNe 

_ N J d£a_ N ' 
y V dw x W 

_ N J dCa N ' 
ly v dv x v 



-N' u T EAeX 

0 

-K T EIyzK 

n: t el z n: +n'jtn' v 


dx 


dx (134) 


0 

0 


NJ ctmA^fNe 
de N e 
0 
0 


-A^csinr^' 


0 

0 

n'J c tan A jJJN' w 

. dc , 


A H T F v GV 0 +A 7 (F*y e + Av* e ) A e 


-A, T csinr|^A„, -A T csinr|^A; 


0 

0 

A^tanA^A,', 

d C m 


dw x w 


-NjFWl 

-Nj {F x y e +F y x e )N' w 

nJfXu, 


-NJ FfN e - A „ 1 (F xZe + F z x e )Ng A „ 1 ( F yZe - F z y e )N w 


> dx 


-NjF y e N v 
Nj (F xZe +F z x e ) A; 
-Al T (F y Ze-F z y e )N' v 

nJfsnJ 


(135) 


(136) 


The effect of the tensile force can be clearly seen in the structural stiffness as the additive stiffness terms in both the 
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bending and torsion stiffness matrices. The presence of a tensile force causes a geometric nonlinear tension stiffening 
effect which increases the natural frequencies of the structural dynamic modes. Conversely, a compressive force causes 
the natural frequencies to decrease. Static instability due to compressive loading could result when the compressive 
force-induced stiffness is large enough to offset the torsion or bending stiffness. 

VI. Applications 


A. Rotary Wings 


Rotary wings are axially loaded structures that experience tensile loading due to the centrifugal acceleration. The 
structural stiffness of a rotary wing is given by 


K s 


N, 


'T 


GJ + EBi (y) 2 

-AC T £fi 2 yA e 

-Nl T EBwN e 


N e -N'jEB 2 yN" v 

< J EI yy K 

-■ K J El yz N'l 


- N e T EB 3 yN ", 

-A CEIyX 

n", t ei zz n". 


dx 


-a 1 


N d J Tk 2 N d 
0 
0 


0 

V T tn' 

iy w ± v w 

0 


0 

0 

n'JtnI 


dx 


where T (x) = J^mrQ?dr is the centrifugal force, r is the radius, and Q. is the angular speed. 

Therefore, the system matrix equations at the angular speed Q. and zero speed can be expressed as 

— ( 0 2 x + (AT 1 K 0 + QrM- l F) x = 0 


(137) 


(138) 


and 


—cofix + M 1 Kox = 0 

Thus, the natural frequencies of a rotary wing are dependent on the angular speed of rotation according to 


(139) 


or = ©a 


(m 2 


(140) 


where fi)„ is the natural frequency at the angular speed Q and (O f) is the zero-speed natural frequency, and / is the 
centrifugal stiffening factor. 

As an example, the natural frequencies of compressor rotor blades in the NASA Ames 1 1 -Ft Unitary Plan 1 1 -Foot 
By 11 -Foot Transonic Wind Tunnel three-stage compressor 9 are computed. The experimental natural frequencies 
were obtained from modal analysis measurements. In addition, wind-on strain gauge measurements were conducted 
at various compressor speeds to obtain the variation in the natural frequencies with the compressor speed. The natural 
frequencies were computed from the Fast Fourier Transform (FFT) of the time signatures of the strain gauges at various 
compressor speeds. This allows the centrifugal stiffening factor to be computed. Table 1 shows the computed natural 
frequencies of the three-stage compressor rotor blades by the finite-element method presented herein as compared to 
the experimental results. The agreement between the computed natural frequencies and the experimental results is 
excellent. 
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Mode 

FEM too (Hz) 

Bench Test ®o (Hz) 

FEM (On at 700 rpm (Hz) 

FEM / 

Wind-On / 

IB 

54.43 

54.18 

61.67 

2.488 

2.375 

2B 

167.32 

160.92 

174.51 

4.252 

4.842 

IT 

302.37 

305.56 

309.25 

5.560 

4.661 

3B 

370.29 

374.61 

375.25 

5.24 

9.110 

2T 

549.98 

554.51 

561.91 

9.87 

8.764 


Table 1 - Natural Frequencies of Rotor Blades in NASA Ames Unitary Plan 11-Foot By 1 1-Foot Transonic Wind 

Tunnel Three-Stage Compressor 


B. Truss-Braced Wing Aircraft 

A recent flutter wind tunnel test of a dynamically scaled TBW was conducted in NASA Langley Transonic Dynamic 
Wind Tunnel (TDT). 3 The model is 16% scaled, semi-span, wall-mounted model as shown in Fig. 3. It was reported 
that the flutter speed is dependent on the angle of attack at ±1 deg. 3 While there may be many nonlinear factors 
that could contribute to this observation, one possible exploration could be the geometric nonlinear tension stiffening 
effect of the main struts at a positive angle of attack, and conversely the softening effect at a negative angle of attack, 
assuming the axial force is below the critical buckling load. 



Fig. 3 - Truss-Braced Wing (TBW) Model in NASA Langley Transonic Dynamic Tunnel (TDT) 

A static aeroelastic analysis is performed to compute the static deflections and the tension carried by the struts. 
The finite-element model (FEM) of the TBW utilizes an equivalent beam/frame configuration that uses a total of five 
components. Only the right side of the aircraft is modeled, and the wing configuration is assumed to be cantilevered 
to the side of the fuselage. The three main physical components of the wing configuration are the wing, the strut, and 
the jury strut. The FEM uses 30 elements for the wing, 28 for the strut, and 5 for the jury strut. Two rigid elements are 
used to connect the strut and the jury strut to the wing. Figure 4 shows the FEM of the TBW. 
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Fig. 4 - Finite-Element Model of Truss-Braced Wing (TBW) 

The FEM uses lumped mass and inertia properties provided by Boeing Research & Technology. 10 The lumped 
mass data is converted into the running mass. Extensional stiffness of the wing is not provided and instead is estimated. 
The engine-pylon mass is not available and therefore is not modeled. This initial model is aimed at understanding the 
geometric nonlinear effect of the tension stiffening, and is not intended to provide an accurate representation of the 
TBW model. Therefore, the results are used only for the purpose of illustrating the geometric nonlinear tension 
stiffness, but are not to be compared to the established TBW data. 

Aerodynamic loads are computed by an aerodynamic analysis code. Figures 5 and 6 are the plots of lift and 
pitching moment distributions for the wing and the strut at the design cruise Cl = 0.81 at Mach 0.7 and an altitude of 
42,000 ft, corresponding to the lg flight load or a cruise gross weight of 170,250 lbs. 


Lift Distribution at 1g 



Pitching Moment at 1g 



Fig. 5 - Lift Distribution Fig. 6 - Pitching Moment Distribution 

A nonlinear static analysis is performed to include the geometric nonlinear effect of tension stiffening. The solution 
is initialized with zero initial tension. After the first iteration, the tension is computed and the stiffness matrix is 
updated. The solution is iterated until the wing tip solution converges. The wing tip deflection from the linear solution 
with no tension stiffening is 1.5434 ft which indicates that the combined wing-strut structure is quite tiff. The nonlinear 
solution predicts a wing tip deflection of 1.5248 ft, which only reduces by about 1%. So the tension stiffening effect 
does not appear to be significant at lg load. For comparison, the nonlinear wing tip deflections at 2.5g and -lg 
are computed to be 3.0622 ft and -1.8885 ft, respectively. These are compared to the corresponding linear wing tip 
deflections of 3.1009 ft and -1.9194 ft. The computed wing static deflections are shown in Fig. 7. 
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Fig. 7 - Wing Static Deflection 


The axial extension of the strut for lg, 2.5g, and -lg flight loads is plotted in Fig. 8. The tension carried by the 
strut is then computed from the axial strain and is plotted in Fig. 9. The static aeroelastic analysis shows that the strut 
on the average carries about 153,000 lbs to support a take-off gross weight of 170,250 lbs. 




Fig. 8 - Axial Extension of Strut Fig. 9 - Tension in Strut 

An eigenvalue analysis is performed to compute the no-load structural dynamic natural frequencies of the TBW 
and the structural dynamic natural frequencies with tension stiffening due to the loaded structure at lg, 2.5g, and -lg 
loads. The natural frequencies of the first ten modes are shown in Table 2. It can be seen that the natural frequencies 
with tension stiffening at lg and 2.5g are larger than the no-load natural frequencies. The increase in the natural 
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frequencies from the no-load natural frequencies ranges from about 1% to 8% for the lg flight load case and from 
about 1% to 15% for the 2.5g flight load case. The -lg flight load case exhibits a significant reduction in the natural 
frequencies ranging from about 1% to 24%, particularly at higher modes. 


Mode 

No-Load 
Frequency (Hz) 

Frequency 
at lg (Hz) 

Frequency 
at 2.5g (Hz) 

Frequency 
at-lg (Hz) 

1 

1.9413 

2.0052 

2.0718 

1.8062 

2 

2.3662 

2.4541 

2.5483 

2.2090 

3 

4.0659 

4.1832 

4.3246 

3.8891 

4 

6.6962 

6.7525 

6.8226 

6.6082 

5 

8.5328 

8.6784 

8.8275 

8.2343 

6 

10.0313 

10.0461 

10.0611 

9.9881 

7 

11.4553 

11.4804 

11.5082 

11.4106 

8 

17.1444 

17.2365 

17.3415 

15.5022 

9 

18.2639 

19.8070 

21.3600 

16.6709 

10 

22.1617 

22.2387 

22.5115 

17.9237 


Table 2 - Natural Frequencies of TBW Aircraft 


A flutter analysis is performed for the no-load case and the three load cases. A 2% structural damping is assumed. 
The no-load case shows a flutter occurring at modes 4 and 5 as these two modes coalesce together, as shown in Fig. 
10. The flutter speeds are 589.00 knots equivalent airspeed (KEAS) for mode 4 and 574.85 KEAS for mode 5, as 
shown in Fig. 11. Figures 12 to 17 are the plots of natural frequencies and damping of the TBW aircraft for the three 
flight loads at lg, 2.5g, and -lg. The flutter speeds for the three flight load cases are shown in Table 3. As can be 
seen, the flutter speeds are dependent on the flight loads due to the geometric nonlinear effect of tension stiffening. 
The difference in the flutter speed is as much as 35 KEAS between the lg flight load and -lg flight load and 15 KEAS 
between the lg flight load and 2.5g flight load. 


No-Load Frequency No-Load Damping 




Fig. 10 - No-Load Frequency 


Fig. 1 1 - No-Load Damping 
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Frequency (Hz) Frequency (Hz) 


Frequency at 1 g Damping at 1 g 




Fig. 142 - Frequency at lg 


Fig. 13 - Damping at lg 


Frequency at 2.5g Damping at 2.5g 




Fig. 14 - Frequency at 2.5g 


Fig. 15 - Damping at 2.5g 
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Frequency at -1g 




Damping at -1g 
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Fig. 16 - Frequency at -lg Fig. 17 - Damping at -lg 


Mode 

No-Load Flutter 
Speed (KEAS) 

Flutter Speed 
at lg (KEAS) 

Flutter Speed 
at 2.5g (KEAS) 

Flutter Speed 
at-lg (KEAS) 

4 

589.00 

601.15 

616.54 

566.48 

5 

574.85 

585.62 

596.83 

555.19 


Table 3 - Flutter Speed of TB W Aircraft 


The flutter Mach numbers at an altitude of 42,000 ft are shown in Table 4. The Boeing Report indicates a dive 
speed of Mach 0.82. A flutter margin of 15% results in a minimum clearance flutter Mach 0.9430. Thus, the flutter 
Mach numbers are well above the clearance flutter Mach. It should be noted that this analysis does not include the 
engine-pylon mass an well as the correct fuel mass. So the computed flutter speeds in this study are optimistic. 


Mode 

No-Load Flutter 

Flutter Mach 

Flutter Mach 

Flutter Mach 

Mach 

at lg 

at 2.5g 

at-lg 

4 

1.2867 

1.3133 

1.3469 

1.2375 

5 

1.2558 

1.2793 

1.3038 

1.2129 


Table 4 - Flutter Mach Number of TBW Aircraft at Altitude of 42,000 ft 


In general, increased wing loading results in an increase in the flutter speed. This is consistent with the observation 
from the wind tunnel test of the TBW in NASA Langley’s Transonic Dynamic Tunnel (TDT) which shows that the 
flutter speed at a negative angle of attack is consistently lower than that at a positive angle of attack at the same 
dynamic pressure. 3 - 10 


VII. Conclusion 

This paper presents a study of the geometric nonlinear effect of axial loading on the aeroelasticity of a truss-braced 
wing aircraft structure. The main struts which restrain the high-aspect ratio wings generally carry a tensile force due 
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to the lift force generated by the wings. The tension causes a geometric nonlinear stiffening effect which effectively 
raises the stiffness of the structure. This geometric nonlinear effect therefore causes the natural frequencies to increase, 
hence changes to flutter speed solutions. A nonlinear finite-element model of the trussed-braced wing aircraft has been 
developed. The nonlinear static deflection analysis is performed and yields a slightly lower wing deflection than the 
linear analysis. Eigenvalue analysis is performed and shows that the natural frequencies at the lg and 2.5g flight loads 
are higher than the no-load natural frequencies. Conversely, the natural frequencies at the -lg flight load is lower than 
the no-load natural frequencies. A flutter analysis is performed and shows that two flutter modes exist. Modes 4 and 
5 flutter at different airspeeds depending on the three flight loads. This is caused by the geometric nonlinear effect of 
tension stiffening due to the strut which affects the flutter frequencies. In general, increased wing loading results in 
an increase in the flutter speed. The study illustrates the importance of accounting for the geometric nonlinear tension 
stiffening effect in analyzing the truss-braced wing aircraft. 
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